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Abstract 

Connectivity patterns of ecological elements are often the core concern of ecologists working at multiple levels of or- 
ganization (e.g., populations, communities, ecosystems, and landscapes) because these patterns often reflect the forces 
shaping the system's development as well as constraining their operation. One reason these patterns of direct connections 
are critical is that they establish the pathways through which elements influence each other indirectly. Here, we tested a 
hypothesized consequence of connectivity in ecosystems: the homogenization of resource distributions in flow networks. 
Specifically, we tested the generality of the systems ecology hypothesis of resource homogenization in 50 empirically 
^ derived trophic ecosystem models representing 35 distinct ecosystems. We applied Network Environ Analysis (NEA) to 
calculate resource homogenization for these models, where homogenization is defined as the ratio of the coefficient of 
variation of the direct flow intensity matrix {CV{G)) to the covariance of the integral flow intensity matrix (CF(N)). A 
ratio greater than unity indicates the presence of homogenization. We also tested the hypotheses that homogenization 
increases with system size, connectance, and cycling. We further evaluated the robustness of our results in two ways. 
First, we verified the close correspondence between the input- and output-oriented homogenization values to ensure 
^T^l that our results were not biased by our decision to focus on the output orientation. Second, we conducted a Monte 
p I Carlo based uncertainty analysis to determine the robustness of our results to ±5% error introduced into the original 
Q flow matrices for each model. Our results show that resource homogenization occurs universally in the 50 ecosystem 
• ^ models tested, with values ranging from 1.04 to 1.97 and a median of 1.61. However, our results do not support the 
hypothesized relationship between network homogenization and system size and connectance, as the results of the linear 
QT^ regressions are insignificant. Further, there is only weak support for the positive relationship between homogenization 
and cycling. We confirm that our results are not biased by using the output-oriented homogenization values instead 
I of the input-oriented values because there is a significant linear regression between the two types of homogenization 
J> (r^ = 0.38, p < 0.001) and the values are weU correlated {S = 8,054, p = 0.61, p < 0.001). Finally, we found that our 
T-H results are robust to ±5% error in the flow matrices. The error in the homogenization values was less than the error 
introduced into the models and ranged from a minimum of 0.24% to a maximum of 1.5% with a median value of 0.58%. 
The error did not change the qualitative interpretation of the homogenization values. In conclusion, we found strong 
. support for the resource homogenization hypothesis in 50 empirically derived ecosystem models. 

Keywords: network environ analysis, indirect effects, input-output analysis, network homogenization, connectivity, 
food web, ecological network analysis 
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"It is a recognized principle of ecology that 
the interactions of organisms and their envi- 
ronment are reciprocal" Redficld (1958) 

1. Introduction 

Connectivity is a core concept in ecology. This is evi- 
dent in Redfield's (1958) opening quote, but other exam- 
ples include Darwin's (1959) entangled bank metaphor, 
exploiter-victim spatial connectivity (Holland and Hast- 
ings, 2008), landscape ecology (Urban and Kcitt, 2001), 
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and community interactions (Jordan et al., 2003). John 
Muir solidified this idea in the public imagination with 
his statement that "When we try to pick out anything 
by itself, we find that it is bound fast by a thousand in- 
visible cords that cannot be broken to everything in the 
universe" (as quoted by Fox, 1981, p. 291). 

Too often this principle idea is dismissively summa- 
rized as: "Everything is connected to everything else" . As 
Peters (1991) pointed out, this statement is vacuous; there 
is no way to determine if this statement is true or false. 
Perhaps more importantly, however, it misses the funda- 
mental point. In ecological systems, as well as other types 
of complex systems, how organisms and their environmen- 
tal components are connected is what is interesting. For 
example, Jordano et al. (2003) suggest that the nestedness 
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pattern commonly found in mutualistic networks may fa- 
cilitate coevolutionary interactions. Dunne et al. (2002) 
found increasing food web connectivity tends to increase 
ecosystem robustness to biodiversity loss, but this effect 
was modulated in part by the pattern of node degree dis- 
tributions. These connectivity patterns reflect the pro- 
cesses that create and constrain a system's development, 
its function and the services it provides. In this paper, we 
investigated a hypothesized consequence of network con- 
nectivity in ecosystems, the systems ecology hypothesis 
that "the action of networks makes the distribution of re- 
sources more uniform" , which is termed network homog- 
enization (Patten ct al., 1990; Fatli and Patten, 1999a; 
J0rgensen et al., 2007). 

The network homogenization hypothesis emerged from 
the development and application of Network Environ Anal- 
ysis (NEA) , one type of ecological network analysis derived 
from the application and extension of Lconticf s (1966) 
economic input-output analysis (see review by Path and 
Patten, 1999b). Patten et al. (1990) noticed that when 
they applied NEA to ecosystem models, the distribution 
of integral flow intensities (the combination of boundary, 
direct, and indirect flow intensities) tended to be more 
even than the distribution of direct flow intensities alone. 
This change in distribution was hypothesized to be a con- 
sequence of the indirect relationships that depend in part 
on the connectivity pattern (Borrctt ct al., 2007, in press.). 

Path and Patten (1999a) introduced a method of quan- 
tifying network homogenization that is now a standard 
part of NEA (Path and Borrctt, 2006; Path and Patten, 
1999b). Subsequently, Path (2004) and Path and KiUian 
(2007) found evidence that network homogenization oc- 
curred regularly in large hypothetical ecosystem networks, 
and it tended to increase with energy-matter recycling as 
well as network size. In contrast to these hypothetical 
models, however, there have been comparatively few ap- 
plications of this analysis to empirically-based ecosystem 
models (but see Path and Patten, 1999a; Borrctt and Os- 
idele, 2007; Gattic et al., 2006; Baird et al., 2009). 

The objective of the work reported here was to to deter- 
mine the evidence for the network homogenization hypoth- 
esis in empirically-based ecosystem models. Specifically, 
we tested two hypotheses. Pirst, we examined the gener- 
ality of the occurrence of network homogenization in 50 
empirically-based trophic ecosystem models. Second, we 
tested the hypothesized relationship between network ho- 
mogenization and network order, connectance degree, and 
the magnitude of recycling. We concluded by considering 
the ecological consequences of this type of connectivity in 
ecosystem networks. 

2. Materials and Methods 

2.1. Models 

For this meta-analysis we used a database of 50 ecosys- 
tem network models that represent 35 distinct primarily 



freshwater and marine ecosystems (Table 1). The mod- 
els exhibit a range of sizes (4 < n < 125), connectance 
(number of direct links divided by the total possible num- 
ber of hnks; 0.03 < (C = i/n^) < 0.40), and recycling 
(0 < FCI < 0.51). At the core of these trophically-based 
networks is a food web; however, any ecological processes 
that transfers energy-matter are considered including nat- 
ural mortality, excretion, and respiration. Purthermore, 
detritus and several other non-living pools of carbon (e.g., 
dissolved organic matter) are usually key components of 
the models. 

These ecosystem models are empirically derived and 
generally meet the network construction criteria of Path 
et al. (2007). We designate them as empirically based 
because empirical estimates were made to quantify some 
number of the energy-matter fluxes; however, they are 
heterogeneous with respect to their degree of empirical 
quantification and resolution. Sometimes to complete the 
model estimates were made from similar systems or gen- 
eralities, rather than the specific system observed (e.g. 
Brylinsky, 1972). Despite this variability, we claim that 
these models are empirically based by contrast to non- 
empirically based ecosystem models such as the hypothet- 
ical models of Webster et al. (1975) or the cyber- models 
built from general ecosystem assembly rules like those of 
Path (2004). 

2.2. Network Environ Analysis and Homogenization 

Our study uses Network Environ Analysis (NEA), which 
is well described in the literature (e.g.. Patten et al., 1976; 
Path and Patten, 1999b; Path and Borrctt, 2006). We in- 
troduce aspects of the output oriented throughflow analy- 
sis most relevant to the work presented. 

NEA is applied to a network model of energy-matter 
storage and flux in ecosystems. In this model, n nodes 
represent species, groups of species, or abiotic components 
and the L weighted directed edges represent the flow of 
energy-matter generated by some ecological process (e.g., 
photosynthesis, consumption, excretion, harvesting, and 
respiration). Let F„xn = (fij) represent the observed flow 
from ecosystem compartment j to compartment i (e.g., 
j — >■ i), Znxi be a vector of node inputs originating from 
outside the system, and yixn be a vector of flows from 
each node that exit the system. This network model is 
like a road map for the transportation of energy or nutri- 
ents through the ecosystem. To apply NEA to the models, 
we usually assume they are at a static, steady-state (bal- 
anced inputs and outputs, but see Pinn, 1980; Shevtsov 
ct al., 2009, for possible ways to relax these assumptions). 
To meet this analytical assumption, we balanced 22/50 
models in our data set that were not initially at steady 
state with the AVG2 algorithm (AUesina and Bondavalli, 
2003). 

Given this model input, our analysis starts with three 
main calculations. Pirst, we determine the throughflow 
vector T, which is the total amount of energy-matter flow- 
ing into or out of each node. This can be calculated from 
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the model information as follows: 

n 
n 



(1) 

(2) 



At steady state, t = (T °"*)^ = T„xi = (T,). From this 
vector, we derive the first whole-system indicator, total 
system throughfiow {TST = X;r=i T). TST indicates the 
magnitude of flow activity in the system and is similar in 
concept to the gross domestic product from economics. 

The second calculation determines the direct flow in- 
tensities, Gnxn — (gij) from node j to i. These are 
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Notice that the elements of G are unitless and that the col- 
umn sums must lie between zero and unity because ecosys- 
tems are open thermodynamic systems (J0rgensen et al., 
1999). 

The final step is to determine the integral flow inten- 
sities, N„xn — (n-ij). These represent the intensity of 
boundary flow that passes from j to i over all pathways of 
all lengths. These values integrate the boundary, direct, 
and indirect flows, and are determined as 
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In equation (4), I = (iy ) = is the matrix multiplicative 
identity, and the elements of G™ are the flow intensities 
from j to i over all pathways of length m. For example, 
the pathway j k ^ i would have a length m = 2. We 
can find the exact values of N because the power series 
converges so that 



N = (I - G) 



(5) 



Multiple manipulations of the G and N matrices have 
been made in NEA (Fath and Borrctt, 2006; Schramski 
et al., 2006), but here we are most concerned with net- 
work homogenization. Fath and Patten (1999a) suggested 
we could use the ratio of the coefficients of variation in 
G and N to quantify the homogenization of resources hy- 
pothesized to occur over the longer pathways captured in 
the integral matrix. Network homogenization is defined as 



HMG 



^ CV{G) 
^ CV{N)' 
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CV{G) and CV(N) represent the standard deviation of 
the elements in the matrices divided by the mean of the 
matrix elements. Fath and Patten (1999a) showed alge- 
braically that Cy(N) scales between and n, but CV{G) 
should be less than n. The HM G ratio has the nice prop- 



erty of indicating how much more evenly the resources 
are spread across the flows when indirect pathways are 
considered relative to the initial distribution in the direct 
flow intensity matrix. If HMG is greater than unity, we 
conclude the network has homogenized the resource dis- 
tribution. We also consider Cy(N) to compare the raw 
differences in resource homogenization between the net- 
work models. 

2.3. Uncertainty Analysis 

To determine the robustness of our results, we tested 
the uncertainty of the NEA homogenization parameter in 
response to potential error in the data used to construct 
the models. We used a Monte Carlo type perturbation 
procedure for this analysis. For each of the 50 models in 
our database, we constructed a set of perturbed models 
by randomly (uniform distribution) altering the values in 
the initial F matrix by ±5%. Following these changes, the 
values of the output vector y were modified as necessary to 
re-balance the models. If a negative y value was required, 
the candidate perturbed model was considered a failure 
because the data must be non-negative for these ecosystem 
models. Candidate models were generated until 10,000 
successful perturbed models were created. We verified that 
the variation in the successful 10,000 models was ±5%, and 
then calculated the NEA homogenization parameter for 
these models. To determine the impact of the ±5% error in 
the flow values, we characterized the resultant distribution 
of the HMG indicator. We chose not to alter the input 
vector z and used a flux perturbation technique that would 
minimize the change to the original model weighted degree 
distributions. 

3. Results 

Our results provide strong evidence for the network ho- 
mogenization hypothesis. The homogenization ratio was 
larger than unity in all models analyzed (Figure la). The 
minimum, median, and maximum values are 1.04, 1.61, 
and 1.97, respectively. This indicates that the elements 
of the integral flow intensity matrix N are relatively more 
similar than the elements of the direct flow intensity ma- 
trix G in all 50 models. Indirect flows distribute the re- 
source flows throughout the system. 

Figure lb reveals that although resource homogeniza- 
tion is happening, a large amount of variability in integral 
flow intensities in the models remains. The coefficient of 
variation in the integral flow matrices ranged from 1.07 in 
the Oyster Reef model to 6.21 in the Lake Quintc model 
after the zebra mussel {Dreissena polymorpha) invasion. 

Figure 2 shows no clear relationship between the de- 
gree of network homogenization and the model number of 
nodes or degree of connectance. Simple linear regression 
for HMG versus model size n {F — 3.65, — 0.07, p = 
0.062) and model connectance C, (F = 0.74, = 0.02, 
p = 0.39) are not statistically signiflcant. Homogenization 
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Figure 1: The degree of resource homogenization in 50 trophically- based ecosystem models. Panel (a) shows the network homogenization 
index HMG, while (b) reports the coefHcient of variation for the elements of the integral flow intensity matrix N, which is the denominator 
of HMG. Error bars in (a) show the range of HMG values for 10,000 perturbations (±5%) of each empirically based ecosystem model. 



does seem to increase with the Finn cychng index FCI; 
however, the shape of the relationship is not weU resolved. 
The relationship does not appear to be linear as a lin- 
ear regression model violates the assumption of residual 
homoscedasticity. FCI and HMG are significantly cor- 
related (Spearman's S ^ 11,355, p = 0.45, p < 0.001), 
but this significance disappears when we exclude the four 
models with an FCI less than 0.01 so it may be driven by 
outlying data points. 

Many recent NEA studies have presented either the 
input- or output-oriented NEA analyses in an effort to sim- 
plify the papers. The authors usually claim the results will 
be qualitatively similar using either orientation. Figure 3 
shows this assumption to be generally true for the network 
homogenization parameter. We found that a statistically 
significant linear regression fits the data (input = 0.69*out- 



put -I- 0.64, p < 0.001). However, it does not explain much 
of the variation as the adjusted was only 0.38. Despite 
this, a Spearman Rank correlation test shows that the val- 
ues are statistically well correlated {S = 8, 054, p — 0.61, 
p < 0.001). 

Our uncertainty analysis suggests the NEA homoge- 
nization parameter is robust to the ±5% variation added 
to the internal system fiuxes F. The error bars in Figure 1 
show the distance between the minimum and maximum of 
the 10,000 HMG values of the perturbed models. These 
ranges tended to be small, with a median and maximum 
value of 0.04 and 0.11. We were unable to generate bal- 
anced perturbed models for the Florida Bay models be- 
cause the output vector y contained 11 initial zero values, 
which consistently led to negative values in the modified y 
and failed perturbed models. 
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Figure 2: The relationship between output oriented network homogenization and (a) the number of model nodes, n, (b) network connectance, 
C, and (c) the Finn Cycling Index, FCI. 
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Figure 3: The relationship between input and output oriented calcu- 
lations of network homogenization in 50 trophically-based ecosystem 
networks. The linear regression is significant with p < 0.001. 



To further estimate the impact of the data uncertainty, 
we divided the inter-quartile range (distance between first 
and third quartiles) by the actual homogenization value 
and multiplied by 100 to get the percent uncertainty. The 
minimum, median and maximum uncertainties for the model 
population was 0.24%, 0.58%, and 1.5% respectively. These 
values are less than the initial ±5% error introduced into 
the F matrix, revealing that the error in the homogeniza- 
tion parameter is generally less than the initial error. De- 
spite the model perturbations, the qualitative interpreta- 
tion of the homogenization parameter does not change. 

4. Discussion 

The core contribution of this work to systems ecology 
and ecological network analysis is that we provide evidence 
for the generality of the network homogenization hypoth- 
esis in empirically derived ecosystem models. This hy- 
pothesis states that fluxes of resources are more evenly 



distributed in the network when integral flow intensities 
(boundary, direct, and indirect) are considered instead of 
just direct flow intensities. This was universally true in 
the 50 empirically derived ecosystem models we examined. 
Further, we found these results to be robust to small per- 
turbations to the original model flux estimates. This sug- 
gests that small uncertainties in the original model data 
would not change the outcome of our analyses. In this 
section, we discuss the context and ecological significance 
of these results and conclude by outlining several possible 
next steps for this research line. 

J^.l. Values of Network Homogenization 

Baird et al. (1991) warn against comparing network in- 
dicators between models with different structures. This is 
problematic because dilferent model aggregation schemes 
as well as other structural features can influence ecological 
network indicators (Calc and OdcU, 1979; Abarca- Arenas 
and Ulanowicz, 2002; AUcsina ct al., 2005; Baird ct al., 
2009). Our goal here is not to compare the individual 
values of HM G, but to provide a broader context to un- 
derstand the values we observed. 

The HMG values for the empirically derived ecosystem 
models (1.04-1.97) tend to be smaller than those found 
in previous studies. Fath (2004) found 2 < HMG < 
2.8 in hypothetical models based on community assem- 
bly rules. This difference in magnitude may be because 
the community assembly models tended to be larger than 
the empirically-based models used in this study, though 
our results do not support the hypothesized relationship 
between HM G and model size. As suggested by Fath and 
Killian (2007), this could also be an artifact of the selected 
pattern of connections in the assembled models. When 
these authors built community assembly type models that 
better matched observed food webs the range of HMG val- 
ues were more similar to those we present. This is further 
evidence for the importance of carefully considering the 
patterns of how organisms and their environmental com- 
ponents are connected. The values previously reported for 
empirically derived models like the Neuse River Estuary, 
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2.07 (Gattie et al, 2006), and Lake Sidney Lanier, 3.10 
±0.31 (Borrett and Osidele, 2007) are also larger than the 
values we found. While there are muhiple factors that may 
cause this including the aggregation decisions, these two 
models are biogeochemically-based networks rather than 
the trophically-based models we analyzed here. This dif- 
ference may be another source of variation (Baird et al., 
2008; Borrett ct al., in press.). 

It is surprising that our results do not provide evidence 
to support the hypothesized positive relationship between 
HMG and model size n and only weak evidence that it 
increases with cycling. We suspect this contradiction oc- 
curs because our model sample is insufficient to show the 
expected pattern. The range of model size remains limited 
(Table 1). However, if this result holds to future scrutiny 
it might suggest an important characteristic of the more 
empirically derived models that is not present in Path's 
(2004) community assembly models. 

The coefficient of variation of the integral flow ma- 
trix, the denominator in the resource homogenization ratio 
(equation 6), shows that despite a tendency for resource 
homogenization, there remains quite a bit of variability in 
the integral flow intensity. This variability increases with 
network size. This pattern is likely driven by the fact that 
the maximum value of CV^(N) increases with model size 
(see Fath, 2004). We don't observe the same trend with 
HMG because CV{G) tends to increase with model size 
at nearly the same rate observed for CF(N). The best fit 
linear regression of CV{G) and n is CV{G) = 0.07*n+1.9 
{p < 0.001, = 0.87) and the linear regression of CV(N) 
and n is CV(N) = 0.05 *n + 0.81 {p < 0.001, = 0.84). 
These trends tend to cancel in the ratio measure. 

4-2. Ecological Significance 

Network homogenization is ecologically significant be- 
cause it concerns the distribution of resources in the sys- 
tem. As stated in the introduction, resource homogeniza- 
tion implies that regardless of where the resources (e.g., 
carbon, nitrogen) enter the ecosystem they are more evenly 
distributed in network models than we might expect from 
the pattern of direct interactions. Presence of network ho- 
mogenization suggests that the indirect flows are distribut- 
ing the model currency (energy-matter) more evenly over 
the links between organisms and their environs. Thus, 
the energy-matter may be more or less available to some 
species than it at first appears. This resource homoge- 
nization may lead to a more distributed control in the 
ecosystems (Schramski et al., 2006, 2007), and transform 
the effective relationships between species such as apparent 
predators effectively operating as a net mutualist (Patten, 
1991; Ulanowicz and Puccia, 1990). 

Resource homogenization is a consequence of the in- 
direct interactions in these networks. Indirect effects are 
hypothesized to generally be the dominant component of 
ecological interactions (Higashi and Patten, 1989; Patten, 
1991; Borrett ct al., 2006, in press.), but this may not 



hold universally. However, network homogenization sug- 
gests that even when indirect effects arc not dominant, 
they may still have important consequences for the eco- 
logical and evolutionary interactions in the system. Fur- 
ther, we wonder if there might be a significant difference in 
HMG between ecosystems classified as healthy and those 
that are stressed or impacted in some way. We hypoth- 
esize that there might be an optimal range of HMG in 
healthy ecosystems. Too little resource homogenization 
and the system vigor or activity is under optimized, but 
too much homogenization and the system becomes crys- 
tallized or brittle in the sense of Rapport et al. (1998) 
and Mageau et al. (1998). This possibility requires further 
investigation. 

4-3. Limitations and Future Work 

Further work in this research line are key to confirming 
the theoretical results. These steps can be divided into 
theoretical and empirical components. 

Theoretical developments are required to address the 
limitations of the work we present. The first issue con- 
cerns our sample of ecological models. While we have 50 
models that represent 35 distinct systems, this remains 
a relatively small sample size that needs to be extended 
as new models become available. Furthermore, many of 
the models were created by a small set of authors. Thus, 
our results could be influenced by their conscious or un- 
conscious modeling biases, including how they choose to 
aggregate species (see Baird ct al., 2009). Given the uni- 
versality of the results, we do not expect this to under- 
mine our conclusions, but it is a potential bias that can be 
addressed in the future. Another limitation of the mod- 
els is that they are all trophically-based models of mostly 
aquatic ecosystems. While we expect the results will hold 
in more biogeochemically-based ecosystem models as well 
as more terrestrial ecosystems, this is a testable hypothe- 
sis. 

Empirical validation of Ecological Network Analysis 
(ENA) like the NEA results reported here generally re- 
mains a challenge for the acceptance and application of 
the theory (Dame and Christian, 2006, 2008). Part of the 
challenge is that many of the ENA predictions are simply 
not directly empirically testable. However, the alignment 
of ENA results like the dominance of indirect effects (Hi- 
gashi and Patten, 1989; Borrett and Osidele, 2007; Borrett 
et al., in press.) and empirical work showing the impor- 
tance of indirect effects (e.g., Menge, 1995; Wootton, 1991, 
1993; Menendez et al., 2007; McCormick, 2009) and initial 
empirical validation attempts such as that by Dame and 
Christian (2008) lend credibility to the work. We specu- 
late, however, that the network homogenization hypothesis 
may be empirically testable using either radio isotopes in 
experiments like Patten and Witcamp (1967) or stable iso- 
topes in a manner similar to Dame and Christian (2008). 
Further, we wonder what correspondence might exist be- 
tween the network homogenization hypothesis and stable 
isotope mixing models. 
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4- 4- Summary 

Our work provides evidence for the generality of the 
systems ecofogy network homogenization hypothesis in em- 
pirically derived ecosystem models. Homogenization was 
universal in the models and the results are robust to po- 
tential data uncertainty. Together with the initial work 
of Fath and Patten (1999a) that introduced the quantifi- 
cation and application of the homogenization metric and 
the work of Fath (2004) and Fath and Killian (2007) that 
investigated the phenomenon in large ecosystem models 
built from community assembly rules, this is strong evi- 
dence for the hypothesized tendency for network organiza- 
tion to homogenize the distribution of resources in ecosys- 
tems. This is an expression of the functional consequences 
of connectivity patterns in ecosystem ecology. 
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Table 1: Fifty empirically derived trophically-based ecosystem models. 



iViOQei 


units 


Tl ' 




1 J 1 




Source 


Lake Findley 


gC m~^ yr^""^ 


4 


0.38 


51 


0.30 


Richey et al. (1978) 


Mirror Lake 


gC m~^ yr~^ 


5 


0.36 


218 


0.32 


Richey et al. (1978) 


Lake Wingra 


gC m-2 yr-l 


5 


0.40 


1,517 


0.40 


Richey et al. (1978) 


Marion Lake 


gC m-2 yr-l 


5 


0.36 


243 


0.31 


Richey et al. (1978) 


Cone Springs 


kcal m-^ yr-i 


5 


0.32 


30,626 


0.09 


Tilly (1968) 


Silver Springs 


kcal yr— 1 


5 


0.28 


29,175 


0.00 


Odum (1957) 


English Channel 


kcal yr— 1 


6 


0.25 


2,280 


0.00 


Brylinsky (1972) 


Oyster Reef 


Kcal yr— 1 


6 


0.33 


84 


0.11 


Dame and Patten (1981) 


Somme Estuary 


mgC d-i 


9 


0.30 


2,035 


0.14 


Rybarczyk et al. (2003) 


Bothnian Bay 


gC yr— 1 


12 


0.22 


130 


0.18 


Sandberg et al. (2000) 


Bothnian Sea 


gC m-2 yr-l 


12 


0.24 


458 


0.27 


Sandberg et al. (2000) 


Ythan Estuary 


gC m~^ yr— 1 


13 


0.23 


4,181 


0.24 


Baird and Milne (1981) 


Baltic Sea 


mgC d-i 


15 


0.17 


1,974 


0.13 


Baird et al. (1991) 


Ems Estuary 


mgC m-^ d-i 


15 


0.19 


1,019 


0.32 


Baird et al. (1991) 


Swarkops Estuary 


mgC d-i 


15 


0.17 


13,996 


0.47 


Baird et al. (1991) 


Southern Benguela Upwelling 


mgC m" ^ d^i 


16 


0.23 


1,774 


0.19 


Baird et al. (1991) 


Peruvian Upwelling 


mgC d^i 


16 


0.22 


33,496 


0.04 


Baird et al. (1991) 


Crystal River (control) 


mgC m-2 d-i 


21 


0.19 


15,063 


0.07 


Ulanowicz (1986); Ulanowciz (1995) 


Crystal River (thermal) 


mgC m— ^ d-i 


21 


0.14 


12,032 


0.09 


Ulanowicz (1986); Ulanowciz (1995) 


Charca de Maspalomas Lagoon 


mgC m— ^ d-i 


21 


0.13 


6,010,331 


0.18 


Almunia et al. (1999) 


Northern Benguela Upwelling 


mgC m" ^ d^i 


24 


0.21 


6,608 


0.05 


Heymans and Baird (2000) 


Neuse Estuary (early summer 1997) 


mgC m^^ d^i 


30 


0.09 


13,826 


0.12 


Baird et al. (2004b) 


Neuse Estuary (late summer 1997) 


mgC m— ^ d-i 


30 


0.11 


13,038 


0.13 


Baird et al. (2004b) 


Neuse Estuary (early summer 1998) 


mgC d-i 


30 


0.09 


14,025 


0.12 


Baird et al. (2004b) 


Neuse Estuary (late summer 1998) 


mgC m— ^ d-i 


30 


0.10 


15,031 


0.11 


Baird et al. (2004b) 


Gulf of Maine 


g WW yr-i 


31 


0.35 


18,382 


0.15 


Link et al. (2008) 


Georges Bank 


g WW yr"! 


31 


0.35 


16,890 


0.18 


Link et al. (2008) 


Middle Atlantic Bight 


g WW m— ^ yr-l 


32 


0.37 


17,917 


0.18 


Link et al. (2008) 


Narragansett Bay 


mgC yr-l 


32 


0.15 


3,917,246 


0.51 


Monaco and Ulanowicz (1997) 


Southern New England Bight 


g WW m— ^ yr-l 


33 


0.03 


17,597 


0.16 


Link et al. (2008) 


Chesapeake Bay 


mgC m" ^ yr-l 


36 


0.09 


3,227,453 


0.19 


Baird and Ulanowicz (1989) 


St. Marks Seagrass, site 1 (Jan) 


mgC m— ^ d-i 


51 


0.08 


1,316 


0.13 


Baird et al. (1998) 


St. Marks Seagrass, site 1 (Feb) 


mgC m d 


51 


0.08 


1,591 


0.11 


Baird et al. (1998) 


St. Marks Seagrass, site 2 (Jan) 


mgC d-i 


51 


0.07 


1,383 


0.09 


Baird et al. (1998) 


St. Marks Seagrass, site 2 (Feb) 


mgC m— ^ d-i 


51 


0.08 


1,921 


0.08 


Baird et al. (1998) 


St. Marks Seagrass, site 3 (Jan) 


mgC m— ^ d-i 


51 


0.05 


12,651 


0.01 


Baird et al. (1998) 


St. Marks Seagrass, site 4 (Feb) 


mgC m-^ d-i 


51 


0.08 


2,865 


0.04 


Baird et al. (1998) 


Sylt R0m0 Bight 


mgC d— 1 


59 


0.08 


1,353,406 


0.09 


Baird et al. (2004a) 


Graminoids (wet) 


gC m-2 yr-l 


66 


0.18 


13,677 


0.02 


Ulanowicz et al. (2000) 


Graminoids (dry) 


gC m ^ yr 1 


66 


0.18 


7,520 


0.04 


Ulanowicz et al. (2000) 


Cypress (wet) 


gC m-2 yr-l 


68 


0.12 


2,572 


0.04 


Ulanowicz et al. (1997) 


Cypress (dry) 


gC m-^ yr-l 


68 


0.12 


1,918 


0.04 


Ulanowicz et al. (1997) 


Lake Oneida (pre-ZM) 


gC m-2 yr-l 


74 


0.22 


1,638 


< 0.01 


Miehls et al. (2009a) 


Lake Quinte (pre-ZM) 


gC m-2 yr-l 


74 


0.21 


1,467 


< 0.01 


Miehls et al. (2009b) 


Lake Oneida (post-ZM) 


gC yr-l 


76 


0.22 


1,365 


< 0.01 


Miehls et al. (2009a) 


Lake Quinte (post-ZM) 


gC m-2 yr-l 


80 


0.21 


1,925 


0.01 


Miehls et al. (2009b) 


Mangroves (wet) 


gC m-2 yr-l 


94 


0.15 


3,272 


0.10 


Ulanowicz et al. (1999) 


Mangroves (dry) 


gC m 2 yr 1 


94 


0.15 


3,266 


0.10 


Ulanowicz ct al. (1999) 


Florida Bay (wet) 


mgC m— 2 yr-l 


125 


0.12 


2,721 


0.14 


Ulanowicz et al. (1998) 


Florida Bay (dry) 


mgC m— 2 yr-l 


125 


0.13 


1,779 


0.08 


Ulanowicz et al. (1998) 



^ 71 is the number of nodes in the network model, C — h j v? is the model eonncctance when L is the number of direct links or energy— matter transfers, 
TST — X] X] fij + X] is the total system throughflow, and FCI is the Finn Cycling Index. 
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